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ABSTRACT 


The  turbomachinery  through-flow  equations  are  reformulated 
for  mass  and  momentum  averaged  quantities.  The  background  of 
this  analysis  is  the  need  for  an  improved-  assessment  of  the 
accuracy  of  through-flow  computations.  Traditional  through- 
flow  analyses  are  based  on  density  weighted  averaged  quanti¬ 
ties  reducing  to  an  area  average  in  incompressible  flows. 

On  the  other  hand,  experimental  data  are  usually  evaluated 
under  the  form  of  mass-averaged  quantities,  particularly  with 
regard  to  the  overall  energy  balance  and  efficiency  estima¬ 
tions.  The  transition  between  these  two  sets  of  quantities 
is  usually  taken  into  account  by  introducing  an  averaged 
aerodynamic  blockage  factor  in  addition  to  the  blade  blockage 
factor  resulting  from  the  density  averaged  quantities. 

The  present  analysis  provides  a  rigorous  derivation  for 
the  momentum  averaged  flow  quantities  and  shows  that  some 
strong  assumptions  on  the  nature  of  the  non-axisymmetric  flow 
components  are  necessary  in  order  to  justify  the  current 
practice  of  introduction  of  aerodynamic  blockage. 

The  recent  availability  of  detailed  flow  data  in  single 
and  two  stage  axial  compressors  allows  a  partial  validation 
of  these  assumptions,  by  the  comparison  of  the  various  non- 
axisymmetric  components.  In  addition  some  guidelines  are 
provided  relating  the  blockage  factors  to  wake  and  loss 


coefficients . 


INTRODUCTION 


Classical  through-flow  analysis  models,  as  applied  in 
aerodynamic  design  systems  of  turbomachinery  are  generally 
based  on  density  weighted  averaged  quantities  (reducing  to 
area  averages  in  incompressible  flows) ,  Smith  (1966) ,  Hirsch 
and  Warzee  (1976),  (1979),  Jennions  and  Stow  (1984) .  This 
formulation  leads  to  the  introduction  of  additional  "inter¬ 
action"  terms,  having  the  same  structure  as  the  turbulent 
Reynolds  stress,  but  arising  from  the  non-axisymmetry  of  the 
flow  and  contributing  to  the  averaged  radial  equilibrium. 

These  interaction  terms  have  to  be  evaluated  explicitly  in 
quasi-three  dimensional  modifications  of  turbomachinery  flows 
whereby  iterative  computations  of  through- flows  (S2  surfaces) 
and  axisymmetric  blade  to  blade  (SI  surfaces)  flows  are  per¬ 
formed,  with  appropriate  data  being  transmitted  from  one  family 
to  the  other  family  of  surfaces,  Hirsch  and  Warzee  (1979), 
Jennions  and  Stow  (1985) . 

An  alternative  to  the  explicit  evaluation  of  the  inter¬ 
action  terms  as  describing  the  effects  of  the  non-axisymmetry 
on  the  averaged  flow,  is  to  introduce  instead,  an  aerodynamic 
blockage  factor,  based  on  mass  flow  considerations,  in  addition 
to  the  blade  blockage  resulting  from  the  density  weighted 
averaging.  This  approach  is  followed  by  Calvert  and  Ginder 
(1985)  in  order  to  define  a  consistent  quasi-3D  interactive 
procedure.  These  authors  rightly  point  out,  that  their  aero¬ 
dynamic  blockage  factor  should  contain  the  same  information 
on  the  non-axisymmetry  of  the  flow  as  the  interaction  terms, 


w  J-_  -P.  ■.-  W  VT^rry , 


and  that  its  use  in  the  continuity  equation  replaces  the  inter¬ 
action  terms  in  the  radial  equilibrium  momentum  equation.  As 
a  consequence  of  this,  Calvert  and  Ginder  define  a  mass-averaged 
through-flow  instead  of  the  density  weighted  area  averaged  flow 
considered  by  the  previously  mentioned  authors. 

Earlier,  both  present  authors  had  separately  stressed  the 
importance  of  the  introduction  of  an . aerodynamic  blockage  in 
through-flow  evaluation  methods,  Hirsch  &  Denton  (1981),  Dring 
(1984).  More  particularly,  this  last  reference  proposes  a 
quantitative  definition  of  blockage  as  the  ratio  between  mass 
averaged  and  area-averaged  axial  velocity  components.  Based 
on  the  extensive  data  base  for  single  and  multistage  axial 
compressors  obtained  in  the  last  years  at  United  Technology 
by  the  second  author  and  his  coworkers,  Dring  et  al.  (1979) , 
(1982),  (1983),  Wagner  et  al.  (1983),  (1984a),  (1984b),  quan¬ 

titative  evaluations  of  the  blockage  factor  were  made  possible, 
showing  its  importance  particularly  near  the  end  walls.  This 
has  been  confirmed  more  recently  by  Dring  and  Joslyn  (1985), 
who  showed  that  both  the  level  and  spanwise  distribution  of 
aerodynamic  blockage  are  important  and  have  a  strong  impact 
on  the  computed  flow  field  at  the  outlet  of  blade  rows.  In 
this  latter  analysis  the  computed  quantities  were  also  con¬ 
sidered  as  mass  averaged  quantities. 

The  debate  between  the  two  families  of  averaged  quantities 
is  central  to  the  validation  of  through-flow  models.  On  one 
hand,  coherent  through-flow  models  can  be  defined  for 


density-area  averaged  quantities,  but  on  the  other  hand 
physical  arguments  and  the  strong  connection  between  mass 
averaged  quantities,  like  stagnation  pressure,  total  energy 
and  machine  efficiency,  are  essential  to  the  correct  estimation 
of  the  energy  exchange  within  the  turbomachinery  blade  row. 

The  present  report  aims  at  the  derivation  of  a  consistent 
through-flow  model  for  mass,  or  more  precisely,  momentum 
averaged  flow  variables  .  A  consistent  model  can  be 
obtained,  at  the  cost  of  six  different  blockage  coeffi¬ 
cients,  depending  on  which  components  of  the  momentum  flux 
are  to  be  considered. 

As  will  be  shown,  if  the  strong  assumption  is  made  of  the 
equality  of  all  the  blockage  coefficients,  then  a  simplified 
model  is  obtained,  which  entirely  justifies  the  semi-intuitive 
approaches  followed  by  Calvert  and  Ginder  (1985)  and  Dring 
and  Joslyn  (1985) . 

Comparison  with  experimental  data  allows  an  evaluation  of 
the  limits  of  validity  of  this  assumption  and  guidelines  are 
presented  for  the  relation  of  the  aerodynamic  blockage  factors 
with  loss  coefficients  in  order  to  enable  its  introduction  in 
design  systems. 

Section  one  will  present  the  recently  derived  averaged 
form  of  the  conservation  equations  in  vector  form,  from  which 
different  formulations  can  be  obtained. 

Section  two  will  discuss  the  important  energy  conservation 
equation  and  the  impact  various  definitions  might  have  on  a 
through-flow  analysis  model. 


The  momentum  averaged  equations  are  derived  in  section 
three  and  the  influence  of  various  blockage  coefficients  is 
investigated. 

1.  DEFINITION  OF  PASSAGE  AVERAGED  FLOW  EQUATIONS 

All  flow  equations  are  averaged  over  the  blade  passage, 
defined  as  the  region  between  the  suction  surface  0  =  0g 
and  the  pressure  surface  of  the  next  blade,  9  =  8  ,  figure 
1.1.  The  area  average  of  an  arbitrary  quantity  is  defined  as 

9 

_  ,  s 

A  =  Q—z  e  /  A  d0  (1. 

p  °s  0 

y  p 

With  the  introduction  of  the  blade  thickness,  d„ in  the 
tangential  direction,  and  the  blade  pitch,  s,  one  can  write 
A  as 

*  ‘  2rrb  J A  dS  < 1  • 

N  P 

where 

»  -  l  -  |  (1. 

and  N  is  the  number  of  blades. 


This  averaging  procedure  is  applied  to  all  the  conservat 
equations  (mass,  momentum  and  energy)  and  in  order  to  handle 


compressibility  effects  a  density  weighted  passage  average 
is  defined,  Hirsch  and  Warzee  (1979) ,  by 


PA 


oA  =  JiSTtJ  J  pA  d6 


(1.4) 


The  deviations  from  axisymmetry  are  defined  by  A’  and  A" 
according  to 


A  =  A  +  A'  =  A  +  A" 


(1.5) 


with 


A'  «  PA"  =  0  (1.6) 

—  (a) 

In  the  following  we  will  use  the  notation  A  instead  of  A, 
ir  order  to  distinguish  those  quantities  from  corresponding 
mass-averaged  values  A^. 


1 . 1  Turbomachinery  Flow  Model 

It  is  important,  in  an  attempt  to  assess  the  validity  of 
different  assumptions  on  through-flow  quantities,  to  keep  in 
mind  the  approximations  at  the  basis  of  the  flow  models 
generally  used  in  turbomachinery. 

The  essential  approximation  is  expressed  by  a 
distributed  loss  model,  Hirsch  (1985)  ,  in  which  the  shear 
stresses  are  replaced  by  a  loss  (entropy)  generating  friction 
force  F^,  considered  as  a  distributed  force,  defined  by  the 


-  ^  ^ 
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total  pressure  loss  coefficients.  Similarly,  the  energy  and 
entropy  equations  are  simplified  by  the  assumption  that  the 
shear  stress  work  is  exactly  balanced  by  the  heat  conduction 
effects,  leading  to  the  following  set  of  equations,  written 
in  the  relative  system 


f|  +  V  (pw)  =  0 


3£-P-W-}  +  V  ( pw  ©  w) 


-  2 p  (oj  *  w)  +  puj  r  +  cF, 


9  (pi) 


+  v  ( pwi; 


'  =  i£ 
at 


(1.7) 


-j^(ps)  +  V(pws) 


w 

T  P  Ff 


In  these  equations,  written  in  conservative  form,  w  is  the 
relative  velocity  related  to  the  absolute  velocity  v  by 


v  =  w  +  u  =  w+toxr 


for  a  steady  rotation  of  angular  velocity 
The  enthalpy , I , is  defined  by 


-2  -2 

T  _  u  4.  W  U 

1  -  h  +  T  T 


=  H  -  U'v 


where  the  stagnation  enthalpy , H , is 


H  =  h  + 


(1.10) 


6 


f f  / 


These  equations  are  considered  as  Reynolds-averaged  for  the 
turbulence  fluctuations,  and  we  will  assume  that  they  are  also 
averaged  for  the  unsteady  effects,  leaving  us  with  a  steady 
flow  in  the  relative  system. 

An  interesting  discussion  of  the  interrelation  between 
different  averaging  scales  in  time  can  be  found  in  Adamczyk 
(1984)  . 


For  steady  relative  flow,  we  have  the  simplified  flow 


model 


V(pw)  =  0 


V  (pw  0  w)  =  -Vp  -  2p(u>  xw)  +  pw  r  +  pFf 


(1.11) 


V(pwl)  =  0 


V  (pws)  =  jjT  p  Ff 


An  alternate  formulation  is  provided  by  Crocco's  form  for 
the  momentum  equation,  where  the  pressure  term  is  eliminated 
by  application  of  the  entropy  relation 


T  ds  = 


dh  - 

p 


(1.12) 


coupled  to  the  non-conservative  forms  for  the  energy  and 
entropy  equations  by  applying  the  continuity  equation.  This 
leads  to 


V(pw)  =  0 


-*•  -*■  -*■  ±  -*• 
-w  x;  =  TVs  -  VI  +  p 


(wV) I  =  0 


T  (wV)  s  =  w  F, 


(1.13! 


It  is  important  to  remember  that  the  entropy  equation  is 
not  independent  of  the  other  equations.  Hence,  when  the 
entropy  equation  is  applied,  as  is  the  case  in  through-flow 
models,  one  of  the  momentum  equations  has  to  be  discarded. 

In  equation  (1.13),  the  absolute  vorticity  X  has  been 
introduced,  defined  by 


C,  =  V  X  v 


1.2  Passage  Averaged  Equations 


(1.14) 


The  above  system  of  equations  are  passage  averaged, 
following  the  definition  (1.2).  A  detailed  derivation  can  be 
found  in  Hirsch  (1984).  As  is  well  known,  the  averaging 
procedure  introduces  a  blockage  factor,  b,  and  a  body  force , 
fB,  as  a  consequence  of  the  three-dimensionality  of  the  flow. 
The  following  system  is  obtained: 


Continuity  equation 


V  ( p  w b)  =  0 


(1.15) 


. \V.S 

•  ,  I  ,  . 

v  v  \ 


f  < 

V  J*  . 

V  o-  fcV 

'Wv 

V>y.' 

I™, 

■w 

•.v. 

■  v  vv 
>>• 


■  V  .* 


>£6 

i  ■**-*•-« 


v.vv 


Moment  tun  equation 


V  ( pw  0  v/b) 


-bVp  + 


-  2 

P  w 


rb  -  2p{u  x  w)b 


+  p?f  b  +  p  ?Bb 


(1.16) 


;--:v 

:v.v 


3S? 

#5 

NT  V* 

%v 

•  v  v 
k 


Energy  equation 


V  ( p  w  I  b  )  =  0 


(1.17) 


The  blockage  factor,  b,  is  defined  by  equation  (1.3)  and  the 
blade  force  is  equal  to 


P  f 


B 


bs  ^Pp 


n 


(P) 


~Ppn 


U) 


)  -  (p 


_  Ee. 


+  p 


s>  7b 
1  b 


(1.18) 


-+  ( p )  (s) 

where  n  c  and  n  are  the  normals  to  the  pressure  and  suction 
surface,  respectively ,  with  components  in  cylindrical  coordinates, 
for  instance,  for  the  pressure  side 


/  -n 


(P)  _ 


39 

_£ 

3r 


—  tan  e ' 
r  p 


-n 


(z)  _ 


39 

3z 


F  tan  Bp 


(1.19) 


(P)  _ 


=  1 


The  lean  angle  and  the  blade  angle  6^  of  the  pressure 
side  are  hereby  introduced  with  similar  relations  defined 
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for  the  suction  side.  The  entropy  equation  will  be  treated 
separately. 


In  cylindrical  coordinates  (r,8,z)  one  obtains  the  following 


projections , 


Continuity  equation 


^(Pwrbr)  +  k(pwzbr)  =  0 


(1.20) 


Momentum  equations 


pwew9 


_  1  i _  ^  A 

r—  -5—  (pw  w  b  r)  +  r—  (pw  w  b  r)  -  - 

br  3r  M  r  r  br  3z  r  z  r 


+  p  ur2  +  2pueu,  *  P(Pfr  +fBr) 


i  fpIP'V'jbr)  4-  i  |j(pwzwzbr)  « 


‘>‘Ff2+£B2> 


(1.21) 


g  _  -t  g  _  M  0 

Br  pwrwgbr)  +  gj-  jglpw^br)  4-  — 


pwawr 


*  -  2c"r“  +  plF£9  +£Be’ 


Energy  equation 


|7(pwrIbr)  +  ^(pw^br) 


(1.22) 


An  interesting  alternative  formulation  for  the  mass  and 


energy  conservation  laws  is  obtained  in  the  axisymmetric 


coordinate  system  (m,n,0),  figure  1-2,  where  n  is  normal  to 
the  meridional  streamline,  m. 

Expressing  mass  conservation  for  the  streamtube  of  thick¬ 
ness  B^,  one  has,  instead  of  equation  (1.20),  since  wn  =  0 
and  assuming  that  the  shape  of  the  meridional  streamline  is 
nearly  axisymmetric, 

-  0  <1-23) 

where  A  =  2-rrrbB^,  and  for  the  energy  equation 

=  0  (1.24) 


In  the  above  expressions,  no  decision  has  been  made  with 
regard  to  the  nature  of  the  averaged  variables  to  be  considered. 
The  current  option  to  be  found  in  the  literature  is  the 
density-weighted,  area  average,  according  to  the  definition 
(1.4),  Hirsch  and  War zee  (1979),  Jennions  and  Stow  (1984). 
Writing  the  density  weighted,  area  averaged  quantities  with 
a  superscript  (a),  instead  of  the  tilda  in  equation  (1.4), 
one  obtains  the  following  consistent  through-flow  model. 


|?(Pbrw^a))  +  f7(Pbrw<a))  =  0 

V(pV(a)  0  w(a)br)  =  -bVp  +  pbF  +  V(x(s)b) 


(1.25) 


7(^(a)I(a)b) 


-  V(pw"l"b) 
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In  equations  (1.25),  the  term  F  is  the  sum  of  all  the  forces. 


F  =  w2r  -  2p(oo  *w^)  +  Ff  + 


(1.26) 


•  53  ( s ) 

and  the  additional  stress  term  t  represents  the  sum  of  the 
"interaction"  terms  expressed  as  the  gradient  of  a  "secondary” 
stress  factor 


=  -  pw"  ©  w" 


(1.27) 


The  energy  equation  shows  that  the  averaged  total  energy. 


I  ^  ,  defined  by 


j  (a) 


h(a)  +  W 


2  (a) 


(1.28) 


is  not  constant  along  a  density  area  averaged  streamline. 

This  has  important  consequences  on  the  consistency  of  through- 
flow  models  and  will  be  discussed  in  more  detail  in  section  2.1. 

1 . 3  Averaged  Crocco's  Form  of  the  Momentum  Equations 

By  averaging  the  entropy  relation  (1.12),  one  can  replace 
the  pressure  gradient  and  the  blade  force,  f0,  by  averaged 
thermodynamic  variables  h^  and  T^  and  by  an  alternate 
blade  force  function  of  the  enthalpy  variations  between  pres¬ 
sure  and  suction  side,  Hirsch  and  Warzee  (1976)  ,  Hirsch  (1984) . 

The  details  of  the  calculations  can  be  found  in  these 
references,  where  the  simplifying  assumption  of  an  axisymmetric 


^  «  .  .  * .  *  .  *%  /  'J  mfj  ^t'7  y  iji  1 1 1,»  iij 


entropy  has  been  made,  in  accordance  with  the  fact  that  the 
friction  force  F^  is  considered,  in  practical  calculations 
as  an  axisymmetric  quantity.  Note  however  that  this  assumption 
can  easily  be  removed  as  seen  in  section  1  of  Hirsch  (1984). 

One  obtains,  for  steady  state  conditions,  with  the  relation 


Vp  -  pf  =  pVh^  -  pT^  Vs  -  p  f 


(i.29: 


where  f,  is  a  body  force  term  expressed  as  a  function  of  h. 


pfh  =  -  ^[(ph"n)p  -  (ph"n)gJ  +  h"Vp  (1.30) 

the  following  expression  for  equation  (1.16) 


j-V  ( pw  ©  wb)  =  p  t  (a)  Vs  -  pVh(a)  +  pFh  (1.31) 


with 


pF 


h 


-  2 

Poo  r 


2p(ojxw)  +  jF.  + 


(1.32) 


When  introduced  in  the  density  weighted  area  averaged 
momentum  equation  (1.25)  one  obtains,  with  the  application 
of  the  continuity  equation 


T (a)  Vs  -  VI(a)  +  Ff  +  ?h  +  -4^(T(S)  b)  (1.33) 

bp 
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In  equation  (1.33),  the  averaged  absolute  vorticity  is  intro¬ 
duced  as 


5<a)  =  5»;la| 


and  the  total  energy  of  the  density  weighted  averaged  flow, I 


(a) 


(a) 


h(a)  +  wU) 
h  +  — =— 


u 


j(a)  _  pv^'w"  H  j(a)  _  £(a)  ( 

2p 

The  quantity  1^  is  the  averaged  total  energy  of  the  flow, 

A  j  v 

while  I  a  is  the  total  energy  of  the  averaged  flow.  These 
two  quantities  differ  by  the  average  of  the  kinetic  energy 
of  the  fluctuations  w" . 


2.  THE  AVERAGED  ENERGY  EQUATION 

The  total  energy  I  (or  H)  plays  an  essential  role  in  all 
the  through-flow  models,  since  the  through-flow  computations 
have  to  rely  on  the  energy  relation  along  the  meridional 
streamlines.  In  most,  if  not  all,  of  the  through-flow  programs 
based  on  axisymmetric  models,  the  constancy  of  the  rothalpy 
is  applied  in  order  to  relate  quantities  in  two  consecutive 
calculation  stations.  When  dealing  with  the  influence  of  the 
non-axisymmetric  effects  it  is  essential  to  be  able  to  esti¬ 
mate  their  influence  on  the  energy  transport  and  exchange. 
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From  the  third  of  equations  (1.25),  one  can  see  that 
neither  the  density  averaged  total  energy  1^,  which 
satisfies 


pij(a>  ti,a) 


m 


3m 


=  -  gV(pw"I"b) 


(2.1) 


(a) 


nor  the  total  energy  of  the  average  flow  I  ,  which  obeys 
the  equation 


pw<*>  i-  i(a> 

m  3m 


=  -  rr9  (  pw"  I"b)  -  £7(pw'Q'k'a;b) 


-  g-V  ( pw"  I  "b)  -  p  w(a)  —  k(a) 
o  m  am 


(2.2) 


are  conserved  along  a  streamline  w^.  On  the  other  hand,  if 

m 

the  streamsurfaces  are  assumed  to  be  axisvmmetric ,  a  unique  mass 
averaged  total  energy  can  be  defined  by 


p«  =  pw(a)I(m) 


1  ,  * 

=  r-  /  (pwl)d9 


(2.3) 


from  equation  (1.14) 


V(p  w(a>  I  (m)  b) 


=  0 


[2.4) 


or  taking  into  account  the  continuity  equation  (1.15) 


pw(a>vT(m)  .  0 


(2.5: 


15 


Inserting  (2.7)  into  (2.6),  one  obtains 


- ,  .  9  2.  (m) 

p  wm  a  35  1 


-  5  w(a>  i-  f("’) 

m  9m 


(2.11) 


showing  that  i'  ,  the  total  energy  of  the  mass  averaged  flow 
is  not  rigorously  conserved  along  a  density  area  averaged 
streamline. 

According  to  Jennions  and  Stow  (1985) ,  the  contribution  of 
the  kinetic  energy  of  the  fluctuations,  can  be  as  high  as  30% 
of  the  total  energy  of  the  averaged  mean  flow. 

Non-axisymmetric  streamsurfaces  -  Radial  mixing  effects 

If  the  blade  to  blade  streamsurfaces  are  not  close  to  an  axi- 
symmetric  shape,  as  a  consequence  of  secondary  flows,  the  assump¬ 
tion  (2.3)  ceases  to  be  strictly  valid. 

Instead,  the  mean  values  appearing  in  equation  (1.22)  are 
evaluated  as  follows,  writing,  see  figure  1.2, 


w  =  w  sino  =  w  sin  (a  +o') 
r  m  m 


(2.12) 


w  =  w  coso  =  w  cos  (o  +o' 
z  m  m 


introducing  hereby  the  averaged  streamsurf ace  slope  angle  a,  and 
the  "twist"  angle  o'  defined  as  the  difference  between  the  actual 
angle  a  and  the  mean  (axisymmetric )  value  o 


o(r,6,z)  =  a(r,z)  +  o'(r,9,z) 


;2.13) 


'•/>* v ».■  «■ '.m  a  .v,*  ^  '.■v  ^  ',Tr??» i 


TT 


Decomposing  (2.12)  and  assuming  that 


cos  a'  =  1  and  sin  o'  =  0 

one  can  estimate  the  averaged  products  as  follows.  For  the  first 
term,  one  would  have 


p  w  I  sin  (a+a ')  =  p  w_I  sino  +  p  w  I  sino'  •  cos  a  (2.14a) 

m  m  m 

and  for  the  second  term 


pwIcos(a+a')  =  pwlcosa-pwlsino'*sino  (2.14b) 

m  mm 

These  expressions  are  introduced  in  the  energy  equation  (1.22). 
The  first  terms  of  (2.14)  give  the  axisymmetric  contribution  of 
the  l.h.s.  of  equation  (2.4).  Indeed,  defining  1^  by 

p  w  I  =  p  w  I  (2.15) 

m  m 


one  has 

fr(P  br  WrI}  +  br  wzX)  =  br  w[Jia)  sin  al(m)) 

+  (p  br  wla)  cos  a  I  (m)  )  +  (b  r  p  w  I  sin  o'  •  cos  c) 
3  z  m  3r  m 

-  I— (b  r  p  w  I  sino'  •  sin  o)  =  0  (2.16) 
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Note  that  the  definition  (2.15)  of  I'  is  identical  to  the 
definition  (2.3)  in  the  case  of  ax i symmetry . 

The  first  two  terms  reduce  to  the  left  hand  side  of  (2.6)  and 
(2.15)  becomes,  with 


p  wm  I  sin  o'  =  (pw^.)  'I'"  =  p*  w^l' 


(2.17: 


where  the  energy  fluctuations  I"'  are  defined  by 


I  =  I(rn)  +  I"' 


(2.13) 


-  -  (a)  3  f(m) 

p  br  w  I 

m  3m 


(b  rp'  w’  I'”  coso) 
dr  r 


(b  rp'w^  I"'  sin  a ) 


b  r  o  '  w  1 


f^br^ 


(2.19) 


where  R  is  the  radius  of  curvature  of  the  average  streamline  m, 
m 

figure  1.2.  The  derivatives  in  the  direction  normal  to  the  axi- 

3 

symmetric  (averaged)  streamsurface  ^  appears  defined  by 


tc  =  cos  a  -  sm  a 
3  n  a  r  o  z 


and  the  curvature  is  defined  by 


The  last  term  of  equation  (2.19)  can  be  neglected  since  the 


gradients  of  the  non-axisymmetric  contributions  are  much  more 


significant  than  their  amplitude.  Therefore,  the  non-axisymmetric 
energy  equation,  generalizing  equation  (2.6)  becomes 


-(a)  3I(m) 

pbrwm  “Si” 


-r— (b  rp'  w'  I' 
3n  r 


(2.20) 


Note  that  the  right  hand  side  of  equation  (2.20)  represents  a 
source  term  originating  from  the  radial  component  of  the  secon¬ 
dary  velocity  field  and  describes  therefore  a  radial  mixing  of 
the  total  energy. 

If  a  gradient  assumption  is  made  for  the  large  scale  non- 
axisymmetric  fluctuations,  one  could  write 


p'  w^I" 

or 


p  '  w'  I"' 
r 

giving  rise  to  a  diflution  type  equation  for  the  energy  redis¬ 
tribution  due  to  the  non-axisymmetric  flow  field. 

When  the  energy  equation  is  written  for  the  total  energy 

a  /  __  \ 

of  the  averaged  flow  I  ,  the  radial  mixing  term  has  to  be  added 
to  the  right  hand  side  of  equation  (2.11),  which  becomes 


3l(m) 

3r 


af(m) 

3n 


(2.21) 


-(a)  3  =(m) 

Wm  aS  1 


—  !m(b  r  o  '  wil"’  ) 
,  —  j  n  r 

brp 


3  =  (m) 

3m  ^ 


(2.22) 
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The  non-constancy  of  the  total  energy  of  the  averaged  flow, 

I  ^  or  1^,  poses  a  fundamental  problem  in  through-flow  compu¬ 
tations  where  the  application  of  some  form  of  total  energy 
conservation  has  to  be  applied.  As  mentioned  above,  the  calcu¬ 
lations  performed  by  Jennions  and  Stow  (1985)  in  the  case  of 
gas  turbine  nozzle  vanes,  indicates  that  the  kinetic  energy  terms 
in  equations  (2.2)  or  (2.11)  are  not  negligible.  One  might  won¬ 
der,  therefore,  what  the  influence  this  might  have  on  the  evalua¬ 
tion  of  the  right-hand  side  terms  in  the  radial  equilibrium 
equation  under  the  form  of  equation  (1.33)  for  instance,  Crocco's 
form. 

The  following  argument  tends  to  support  the  statement  that 
the  entropy  variations,  or  more  precisely  the  rotary  stagnation 
pressure  gradients,  are  more  important  than  the  enthalpy  varia¬ 
tions,  at  least  for  low  speed  flows. 

From  the  isentropic  relations  between  static  and  absolute 
or  relative  stagnation  conditions,  one  has 

T  ds  =  dh  -  ^ 

p 


or 


T0  ds 


where  the  subscript  indicates  absolute  stagnation  conditions. 
For  rotary,  relative  conditions,  whereby 
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T  + 


(2.23) 


Hence 


T* 

0 


2c 

P 


2c 


one  has 


p(^A-l 


ds 


=  dl  - 


dp 


Tds  -  dl  » 


dl - —  dp*  -  dl 


*_  * 
P0T0 


(~T  ~  1)  dl 


El  dp* 
p 

0 


<-7-l)dI 

T0 


pPr 


7  dP, 


(2.24) 


Introducing  a  rotary  Mach  number,  M 


0 


M 


★  2 


-*■2  -*-2 
w  -  u 

yrT 


(2.25) 


the  coefficients  in  equation  (2.24)  become, 
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The  terms  (TVs  -  VI)  in  the  right-hand  side  of  the  momentum 
equation  (1.33),  becomes,  independently  of  the  type  of  averaged 
quantity  considered 


TVs  -  VI  = 


i  +SEI 

2 


(1  +1i1m*2)y'/y"1  ^P° 


(2.26) 


In  low  speed  compressors  the  blade  exit  cotary  Mach  number 
(2.25)  should  be  small.  This  is  in  particular  the  case  for  the 
data  of  the  UTRC  experimental  axial  compressor  runs  discussed  in 
the  next  section.  When  this  is  the  case,  the  influence  of  the 
stagnation  pressure  variations  is  clearly  dominating  the  through- 
flow  and  the  momentum  exchange. 

One  could  therefore  consider  that  small  errors  on  the 
evaluation  of  the  enthalpy  ’  transport  will  not  affect  signifi¬ 
cantly  the  radial  equilibrium  of  the  flow. 
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3.  MOMENTUM  AND  MASS  AVERAGED  THROUGH-FLOW  EQUATIONS 

With  density  weighted,  area  averaged  variables  one  can 
clearly  define  a  consistent  through-flow  model  where,  next 
to  the  blade  blockage  factor  b,  interaction  terms  due 
to  the  non-axisymmetric  secondary  stresses  x^  are  describing 
the  three-dimensional  effects  on  the  average  flow. 

As  mentioned  in  the  introduction,  mass  averaged  variables 
such  as  stagnation  pressure  and  total  enthalpy  are  more  repre¬ 
sentative  of  the  physical  energy  exchange  than  the  corresponding 
area  averaged  variables.  In  an  attempt  to  formulate  a  coher¬ 
ent  through-flow  model  for  mass-averaged  variables  we  recon¬ 
sider  the  passage  averaged  equations  (1.20)  to  (1.22),  where 
for  simplicity  the  momentum  equations  are  considered  in  the 
absolute  frame  of  reference. 

The  mass  conservation  equation  and  the  energy  equation 
have  already  been  discussed,  with  the  following  outcome. 

Mass  conservation 

The  natural  averaged  quantities  are  the  density  weighted, 
area  averaged  velocities  which  lead  to  the  first  equation 
(1.25),  reproduced  here 

!^(pbrw^)  +  |^(pbrw^a))  =  0  (3.1) 

or  in  vector  form 

V(pbrwU))  =  0  (3.2) 
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3.1  Momentum  Equations — Momentum  Averaged  Velocity  Components 


Considering  equations  (1.21)  in  the  absolute  system,  it  appears 
that  it  is  not  possible  to  define  a  unique  momentum  averaged  ve¬ 
locity  component,  since  different  momentum  flux  components  appear 
in  the  projections  of  the  equation  of  motion.  For  instance,  in 
the  radial  component,  the.  averaged  flux  components  pvrv  ,  ^vrwz 
and  pvTvT  occur,  while  in  the  axial  projection  one  encounters  the 

U  U 

components  Pvzvr  and  pvzvz .  This  leads  to'  the  definition  of 
momentum  averaged  velocity  components  such  as 


pv  v 

r  r 


~-v(r) 
r  r 


p  ?(r,7(a) 
r  r 


(3.3) 


pv  v 
r  r 


E  f  0Vr  Vrd8 
P 


(3.4) 


pv  v 
r  z 


—  ?,r) 

r  z 


5  v<a) 


(3.5) 


pv  v 

r  z 


,  s 

b  )  p  vr  vz  d0 
P 


(3.6) 


One  has  also 


pv  v 

r  z 


- (a)  =(z) 

p  vz  vr 


—  (  ) 

The  averaged  velocity  component  v^  represents  the  average  of  v^. 

weighted  by  the  radial  momentum  pvr-  Similarly,  from  (3.6)  and 

(3.7),  one  can  define  another  averaged  radial  velocity  component 
—  (  2  )  4 

vr  ,  representing  an  average  vr,  weighted  this  time  by  the  axial 
momentum  pv  .  A  priori  there  is  no  reason  to  consider  these  two 


components  as  equal.  Prom  the  definition  of  the  velocity  fluctua¬ 


tions  v",  following  equation  (1.5),  one  can  write  for  instance 
for  pvrvz. 


pVz 


p  v<a)  V{a) 
r  z 


+  pv”v" 
r  z 


(3.8) 


where  pv"v"  are  the  components  of  the  secondary  stress  tensor 

IT  2 

=  ( s ) 

t  considered  in  the  absolute  system. 

We  define  now  six  different  blockage  coefficients,  K^_.  where 
i  and  j  represent  the  components  (r,9,  z)  by 


pvivj_ 


—  — (a)  —(a) 
p  v  .  v  . 
i  3 


_J _ 

-(a) 
v  : 

3 


(3.9) 


Hence  K~^  is  the  ratio  between  the  averaged  velocity  component 
Vj  weighted  with  the  momentum  pv^  divided  by  the  corresponding 
density  weighted,  area  averaged  velocity  component  v.ja^ • 

From  the  symmetry  between  i  and  j  best  seen  in  equation 
(3.10)  one  has  also 


i  “h1  i 

K .  ^  =  -7— r  *  K.  1 
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(3.10) 


or 


K.  . 
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=  (j) 


v.(a) 

1 


(3.11) 


This  definition  of  a  blockage  coefficient  with  respect  to  a  selected 
momentum  component  generalizes  the  definition  of  the  aerodynamic 
blockage  coefficient  introduced  in  Dring  (1984)  and  discussed  more 
in  details  in  Dring  and  Joslyn  (1985).  In  these  references,  the 


is  determined  from  the  mass  averaged 


averaged  velocity  v^ 


dynamic  pressure,  and  mass  averaged  flow  angles  assuming  constant 
flow  angles  equal  to  their  averaged  values. 

From  equation  (3.9)  one  can  also  write 


pvVv'.’ 


K-1  =  1  + - (a)  ^  —  ( a ) 

11  p  v .  v : 
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(3.12) 


(7^3) 


Introducing  the  v:J  component  in  the  equations  (1.21)  leads  to 
the  following  formulation  for  the  radial  component 


bf  !?(p  vra)  vir)br)  +  '^rZ)br)  - 


4*  +  P(Ffr  +  fBr> 


3r 


(3.13) 


or  in  function  of  the  ^  coefficients,  eliminating  the  density, 


area  averaged  component. 


55  Krr  v^’v^’br)  +  55  Krz 


o  K  ^(r)  v{9) 
p  K9r  9  v9 


'  If  +  p(Ffr  +fBr) 


:  3 . 14 ) 


An  alternative  formulation  is  obtained  as  a  function  of  the 


density  weighted  velocity  components  v^a^ .  Equation  (3.13! 


becomes  by  eliminating  the  variables  vf^. 
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**.  ■*»  «*/**J*^ . v  * /  v  v  v 


V  \  s'  %  v 


'  ^  ^  %  v  V  V  V 


T~-  I — (p  V(a)  v(a)brK_1)  +  t—  |-(p  v(a)  v(a)brK_1) 
br  3r  r  r  rr;  br  3z  z  r  rz 


p  H*’  ^  _-i 

- r - Kee 


3r  +  p(Ffr  +  fBR} 


(3.15) 


Equation  (3.15)  is  to  be  compared  to  the  projections  of 
equation  (1.15);  for  instance  the  radial  component  gives  in 
the  absolute  system 


£fe(^>v<*>br)  |j(p7U>v<*>be)  - 


+  P(FC  +  f_  )  -  xr-  r(Pv'V'br) 
3r  fr  Br'  br  3r  r  r 


.  ov"v" 

-  si  3?‘^br>  *  -hr 


(3.16) 


It  can  be  seen  from  the  relation  (3.12)  that  both  formulations 

are  identical  and  that  the  coefficients  do  contain  the 

•  »  =  ( s ) 

same  information  as  the  secondary  stress  x  .  Hence  the 

interaction  terms  do  not  appear  in  equation  (3.15),  since  the 

whole  influence  of  the  non - ax i symmetry  is  contained  in  the 

K—  coefficients.  The  other  components  of  the  momentum 

equation  can  be  treated  in  a  similar  way. 


Axial  momentum  equation 

One  obtains,  from  the  second  equation  (1.21), 


~  f-(p  v{a)  v[z)br)  +  ~  |—  ( p  v(a)  v(z)br) 
br  3r  z  r  br  dz  z  z 


If  +  p(Ffz  +  f B; 


(3.17) 
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The  alternative  forms  as  a  function  of  the  momentum  and  density 
averaged  velocity  components  are  as  follows. 


E7  vzr)  vrr>Krrbr>  *  SF  !?(p  vz2’Krzbr) 


"  '  H  +  p(Ffz  +  £Bz> 


:3. is) 


i  3  (p  =(D  7‘2»k  br).  +  pir  |r(s  7<2)  V ' z > K  br) 
br  dr  z  r  zz  br  dz  z  z  zz 


=  ~  If  +  P(FfZ  +  fBZ) 


3.19) 


In  function  of  the  v^a^  components  one  has 


br  dr 


v(a)  v^a)K  ^br)  +  -L  i_(p  v(a>  v (a)  K~^br) 
z  r  rz  br  3z  z  z  zz 


=  “  If  +  p(Ffz  +  fBz} 


(3.20) 


where  again,  by  comparison  with  equation  (1.25)  the  interaction 

terms  are  fully  absorbed  by  the  K  \  and  K  ^  coefficients. 

r  z  z  z 


Tangential  momentum  equation 


,v<a)v*r) 

9  o 


p(Ff0+  f B0 ' 


3.21) 


y-y- 


y.-y-l 


.NV-V'.'-V-1 


s'  s'  s'  s'  s'  s' 


*  4"11  ■*-  m  j  t .  »  _  V  *. 


Replacing  the  components,  leads  to  the  following  forms. 


^f?(p-rr)  v9r)Krr*r)  -  ^|l^^r)veZ)Krzbr) 


P^(6)^r) 

+  __  -  KrQ  =  P(FfQ+  f  BQ ) 


(3.22). 
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~1 


or  introducing  other  variables 


Er  f?(p  vr  ^  v0Z)Kz0bl:)  +  EF  ll(s  vzSI  v(Z)Kz9br) 


p  =<9>  7«b> 

+  r  Kr0 


P(Ff0  +  fB9) 


(3.23) 


'  N"  I 
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In  function  of  the  density  averaged  velocities  one  has 


br  f?(p  vra>  v6a)  Kr0brl  +  E7  fe‘p  ^  v0a,K'z0br) 


>V>1 


»  -  1  Jr* 


.  p  E.i 

r  Kr0 


p  (Ff 9  +  fB0) 


(3.24) 


to  be  compared  with  the  third  of  the  momentum  equations  (1.25). 


Continuity  equation 


^•->2 


In  a  similar  way,  the  continuity  equation  can  be  rewritten 
as  a  function  of  the  momentum  averaged  velocity  components,  in 

different  ways,  according  to  the  choice  of  the  coefficients. 

From  equations  (3.1)  and  (3.11)  one  can  write,  with  the  absolute 


m 


30 


velocity  component 


3,^ 


K  br) 
rr 


3 

3l(p  V 


~=(r) 
2 


K  br) 
zr 


=  0 


(3.25) 


or 


37(p  vr 


K  br) 
rz 


4-(p  v!z| 
d  Z  2 


K  br) 
zz 


(3.26) 


A  difficulty  appears  when  the  equations  have  to  be 
derived  in  other  forms,  for  instance  in  the  non-conservative 
form  or  Crocco's  form.  Since  one  cannot  introduce  the  con¬ 
tinuity  equation  in  the  conservative  form  of  the  averaged 
momentum  equation  as  a  consequence  of  the  non-equality  of 
velocity  components  such  as  v^r^  f  v^2^ . 

Indeed,  the  left  hand  side  of  equation  (3.13)  can  be 
worked  out,  leading  to 


p7U)  =<r> 

w  r  dr  r 


- (a)  3  — ( z ) 

p  v  v 

z  3z  r 


pr(v(r) 

br  L  r 


3  , - (a) 

vr 


br) 


+ 


br)  ] 


r 


fBr} 


(3.27) 


If  v^r)  =  v^2)  one  could  factor  out  this  velocity  component, 
and  the  term  in  brackets  would  vanish  due  to  mass  conservation. 
Otherwise,  one  could  not  obtain  Crocco's  form  in  a  consistent 
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way,  without  the  simultaneous  appearance  of  different  com- 

w  =(z)  =(r) 

ponents  vr  '  vr 

3 . 2  Simplifying  Assumptions--Unigue  Blockage  Coefficients 
The  assumption  is  made  of  a  unique  ratio  between  density 
and  momentum  averaged  quantities.  This  corresponds  to  the 
assumption  that 
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(3.28) 


implying  that  all  the  coefficients  are  equal.  Hence,  if 

this  is  satisfied,  there  is  a  unique  definition  of  mass- 
averaged  velocities,  since  for  the  axial  component  for 
instance , 


I _ J 


:::sj 


v  (m) 

z 


pvzVz 

-(a) 

pvz 


pv  v 
^  z  m 

—  (a) 


=  K'1  v(a) 
z 


(3.29) 


and  similarly  for  the  other  components. 

With  this  assumption,  the  continuity  equation  (3.25) 
or  (3.26)  becomes 


o-.'.-N 


|— (p  V(m)  Kbr )  +  |-(p  vjm)Kbr)  =  0 


(3.30) 
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This  equation  is  similar  to  the  density  weighted  formulation 
(1.25)  where  the  replacement  by  mass-averaged  velocities 
has  led  to  the  introduction  of  a  global  blockage  coefficient, 

B  =  Kb.  This  coefficient  is  a  product  of  the  total  blade 
blockage  and  of  an  aerodynamic  blockage  K  due  to  the  no n- 
axisymmetry  of  the  flow. 

The  momentum  equations  can  be  reduced  as  follows,  consider¬ 
ing  first  the  radial  component  (3.14). 

With  the  continuity  equation  (3.30),  one  obtains 


PK[v^> 


-<m> 
dr  r 
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3_ 

3z 


-<m) 

r 


J 


pv 


=  (m)  =  (m) 


r 


K 


fBr> 


(3.31) 


Introducing  equation  (1.29)  in  the  right  hand  side  and  writing 


r 


9  ,  — (m) , 
5r lrve  1 


one  obtains,  with  the  assumptions  that  the  thermodynamic 
variables  T^  and  h^  relate  to  their  mass  averaged  values 


in  the  same  way  as  the  velocity  components. 


=  (m) 

=  (m)  3_  =(m)  +  =(m)  3_  =  (m)  _  _J_  i_(rv<m>)  +  v(m)  i-  v(m) 

r  3r  r  z  3z  r  r  3r  9  9  3r  9 


T(a)  3s  1  3h (a)  Ffr  ffr 

K  or  ~  K  3r  K  K 


f(m)  |2.  -  I-  Eiah  *  Ff-r-  ■ +  h(a)  Jsli 

3r  3r  '  K  or 


=  (m)  3s  3 _ ff  ( m)  ,  (v(m))  ,,13  ,-(m)v2 

T  37  '  37th  +  - 2 - ]  +  2  37(v  5 


Ffr  +  f Br  ,  c(a)  3k"1 
K  3r 


(3.32) 


Finally  one  obtains 


-(m)  r  3  =(m)  3  -(m) 

V  (T —  v  -  - —  V„ 

z  3z  r  or  z 


=  (m) 

9  3  ,  =(m) . 

—  a7(rve  > 


_  =(m)  _  _9_  h(it^  +  F£r  +  fBr 

3r  3r  K 


h(m)  ^  in  K  (3.33) 
or 


Within  the  above  simplifying  assumptions,  the  through-flow 
equations  can  be  interpreted  as  functions  of  the  mass-averaged 
flow  variables.  The  influence  of  the  "interaction"  terms  in 
equation  (1.25),  is  taken  up  by  the  K-coef f icients . 

In  the  continuity  equation  the  geometrical  blockage  factor 
b  =  1  -  d/s  is  to  be  multiplied  by  K,  leading  to  an  overall 
blockage  coefficient 
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B  =  bK 


which  has  to  be  taken  into  account  also  in  the  axial  regions 

between  adjacent  blade  rows,  e.g.,  where  b  =  1,  but  K  jt  1. 

The  additional  term  [h  ^  I—  in  K]  can  be  added  to  the 

d  r 

force  terms  [Ffr  +  fgr3  divided  by  K  and  form  a  generalized 
force  component.  Remember  that  f^  also  contains  a  contribution 

_  g 

equal  to  p  gY  in  b. 

Equations  (3. 30)  and  (3.33)  are  consistent  with  the  radial 
equilibrium  treatments  of  Calvert  and  Ginder  (1985)  as  well  as 
Dring  and  Joslyn  (1985)  ,  in  terms  of  mass-averaged  quantities 
and  a  simple  aerodynamic  blockage  coefficient  replacing  the 
interaction  terms  in  the  density  weighted  formulation. 

The  energy  equation  follows  from  (2.22)  or 


=(m)  3  l(m) 
w  —  I 
in  om 


wrim)  k(m)  +  —  Ijrtbrp'w'I1"  ) 

m  3m  pbr  3h  r 


(3.34) 
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In  a  non-rotating  blade  row,  the  above  equation  reduces  to 


=  (m)  3 _  ~(m) 

m  3m  H 


w'm)  k(m)  +  —  |r-(brp’w;H’,'  ) 
m  dm  — ,  on  r 

pbr 


(3.35) 


showing  that  Hu  as  appearing  in  equation  (3.33)  is  not  con¬ 
served,  even  in  a  non-rotating  blade  row.  However  as  discussed 
in  section  2,  in  combination  with  the  entropy  term,  the 
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dominating  effects  will  most  probably  come  from  the  stagnation 
pressure  variations. 

Note  that  the  tangential  equation  (3.22)  becomes 


-  =(m) 


m  3  r  =(m) , 
—  3HIrvo  1 


°K<Ff0  +  £B6 1 


(3.36) 


Comparison  With  Data 

From  the  data  base  collected  in  the  last  years  at  UTRC, 
the  coefficients  can  be  evaluated  at  different  stations 
of  a  single  and  two-stage  axial  compressors. 

Figures  3.1  to  3.3  show  the  difference  between  mass-averaged 
and  area-averaged  static  and  absolute  and  relative  stagnation 
pressures  at  the  exits  of  stator  1,  rotor  2  and  stator  2  of 
the  two-stage  compressor  reported  in  Dring  et  al.  (1982),  (1983). 

The  common  observation  is  the  nearly  axisymmetric  behavior 
of  the  static  pressure,  while  the  stagnation  pressures  show 
large  non-axisymmetric  effects  in  the  end-wall  regions. 

Figures  3.4  to  3.6  show  the  corresponding  ratios  of  mass 

to  area  averaged  absolute  or  relative  velocities  squared,  as 

well  as  the  ratios  of  the  averaged  kinetic  energies.  The 

curves  follow  very  closely  the  stagnation  pressure  variations, 

as  expected  since  static  pressure  is  nearly  constant.  In 

addition,  the  difference  between  the  two  curves  is  a  measure 

of  the  influence  of  the  kinetic  energy  of  the  fluctuations, 

k^  .  This  influence  is  mostly  sensible  in  the  end  wall  regions 

2 

as  can  be  seen  from  figure  3.7  where  the  quantity  [1  +k 
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is  plotted  in  function  of  space.  Figures  3.8  to  3.10  display 

the  K. .  coefficients  at  the  same  three  locations. 
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The  validity  of  the  uniformity  assumption  (3.28)  can  be 
estimated  on  the  basis  of  these  results.  For  the  components 
not  involving  the  radial  velocity,  the  assumption  of  equality  of 
the  is  satisfied  with  an  acceptable  accuracy  even  in  the 

end  wall  regions. 

The  results  of  the  coefficients  involving  the  radial  velocity 
components  are  more  puzzling.  It  should  be  noticed  however  that 
due  to  their  small  magnitude,  a  large  error  is  connected  to  the 
determination  of  the  Kr^  coefficients,  and  more  data  would  be 
required  from  turbomachines  having  larger  average  radial  velocities. 


4.  ESTIMATION  OF  BLOCKAGE  FACTOR 

If  the  blade  to  blade  flow  can  be  separated  into  a  wake 
and  an  inviscid  region,  one  can  easily  derive  a  relation  between 
the  unique  blockage  factor  and  boundary  layer  or  wake  parameters, 
implying  a  connection  with  total  pressure  loss  coefficients. 

In  figure  4.1  from  Dring  and  Joslyn  (1985),  the  above  assump¬ 
tion  is  valid  for  the  flow  at  exit  of  the  rotor  at  all  span 
locations,  with  the  exception  of  the  tip  clearance  flow  which 
completely  perturbs  the  blade  to  blade  distribution.  This  is 
confirmed  by  the  single  and  2-stage  data. 

If  02  is  t-he  outlet  flow  angle,  the  boundary  layer  thickness 
is  measured  in  the  direction  normal  to  the  flow,  that  is  a  dis¬ 
placement  thickness  5*  is  defined  by,  see  figure  4.2- 


/  (i  — £¥— )  dn 

pJ  peve 
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p  e  e 


(4.1) 
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bs  cos  B- 


P  v 
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bs  cos  B- 


Note  that  this  thickness  5  contains  the  contribution  both  from 

the  section  and  pressure  surfaces,  <5*,  5*  and  the  trailing  edge 

p  s 

thickness  tTE* 

Similarly,  a  momentum  thickness  0,  is  defined  by 


s  cos  3- 
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P  V 
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The  density  averaged  absolute  velocity  is  related  to  5  by 
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p  v 
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7h  /  p  vdv  = 


pe  ve  ^ 1  bs  cos  S2  ^ 


(4.4) 


If  oe  is  taken  equal  to  the  averaged  density 


pe  "  p 


(4.5) 


and  defining  a  dimensionless  displacement  thickness  A^  by 


bs  cos  B. 


^+6s+tTE 
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one  has 


V(a)  =  v  (s  -  A, ) 
e  l 


(4.7) 


The  mass-averaged  velocity,  can  be  estimated  as  follows 
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or,  assuming  a  uniform  flow  angle  in  agreement  with  the 
hypothesis  of  a  unique  blockage  coefficient  K, 


to  the  total  pressure  loss  coefficient  - _  where  q 

q 

is  an  appropriate  dynamic  head.  When  q  is  taken  as  the 
downstream  dynamic  head  in  the  mixed-out  flow,  then 


(4.14) 


and 


/l  -  (Apt/q) 


(4.15) 


which  is  the  relation  derived  in  Dring  (1984)  . 

These  relations  should  allow  the  designer  to  estimate 
the  blockage  factor  in  relation  with  the  assumed  loss 
distribution. 


Conclusions 

A  detailed  analysis  of  the  momentum  averaged  formulation 
has  shown  that  a  simplified  model  can  be  obtained,  in  the 
line  of  more  intuitive  considerations,  if  a  simple  blockage 
coefficient  can  be  assumed. 

This  coefficient  is  then  defined  as  the  ratio  between  the 
density  weighted,  area  averaged  velocity  and  the  mass 
averaged  velocity.  It  contains  the  full  information  of  the 
influence  of  the  non-axisymmetry  of  the  real  flow  on  its 
averaged  through-flow. 
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The  present  analysis  therefore  confirms  the  importance 
of  the  aerodynamic  blockage  as  a  major  parameter  to  be  intro¬ 
duced  to  the  through-flow  design  systems,  instead  of  the 
interaction  terms  based  on  the  secondary  stresses. 

When  the  flow  can  be  separated  into  a  wake  region  and  an 
inviscid  zone,  a  relation- can  be  established  between  this 
blockage  coefficient  and  wake  parameters  or  loss  coefficients 
This  is  however  not  possible  when  strong  leakage  and  strong 
end-wall  stall  is  present  and  no  wake  regions  can  be  defined. 
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Figure  3.7  -  Comparison  of  ratio  of  mass-averaged  and  area 
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